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We develop a dipole model for HERA DIS data which incorporates the impact parameter dis- 
tribution of the proton. The model describes the inclusive total 7*p cross-sections as well as the 
diffractive J/tp differential cross-sections. We compare the model with previous approaches and 
show that the t-distributions are sensitive to saturation phenomena. We estimate the boundary of 
the saturation region and show that it dominates the data in the low-Q^ region where the total 
7*p cross-section exhibits the same universal rise as hadronic cross-sections. The model is then 
extended to nuclei and shows good agreement with the nuclear shadowing data at small- a:: . Finally, 
we estimate the saturation scale in nuclei. 
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One of the most important observation of the HERA 
experiments is the rapid rise of the 'y*p cross-section with 
the 7*p center of mass energy W in the deep inelastic 
scattering (DIS) region. This observation indicates that 
over a region of size ~ l/Q^ additional gluons are abun- 
dantly radiated at small x w /W^. The simultaneous 
observation of a significant number of diffractive events in 
DIS lead to the dipole saturation model of Golec-Biernat 
and Wiisthoff (GBW) The model successfully de- 
scribed the data in the transition region between low and 
intermediate and indicated the emergence of satura- 
tion phenomena in DIS. In spite of its compelling success 
the model has some obvious shortcomings. First, the 
treatment of QCD evolution is only rudimentary. Second, 
the dipole cross-section is integrated over the transverse 
coordinate although the gluon density is expected to be a 
strongly varying function of the impact parameter. The 
evolution was improved by incorporating DGLAP into 
the dipole cross-section thereby extending the model ap- 
plicability to the high region . 

In this paper we further improve the dipole model by 
introducing the impact parameter of the proton into the 
dipole dynamics. For the form of the dipole cross-section 
we adopt the Glauber-Mueller formula which can be de- 
rived under simplifying assumptions 3j . The transverse 
profile of the proton can be probed experimentally by 
measuring the momentum transfer t to the proton in 
diffractive events. Thus diffractive and inclusive j*p pro- 
cesses form a tool to determine the gluon density and sat- 
uration dynamics at HERA. Whether the experimental 
data are sufficient to reveal the rich non-linear dynamics 
of saturation and the Color Glass Condensate jjilremains 
an important phenomenological question HQl which 
this work addresses. 

The Impact Parameter (IP) Saturated dipole model 
can be extended to DIS interactions in nuclei in a nat- 
ural way. We study nuclear effects and compare them 




FIG. 1: 

proton. 



The interaction of a quark anti-quark dipole with a 



with data on nuclear shadowing. Finally we estimate the 
saturation scale in heavy nuclei. 



II. DIPOLE DESCRIPTION OF DIS 

In the dipole picture 0|, the j*p interaction pro- 
ceeds in three stages: First, the incoming virtual pho- 
ton fluctuates into a quark-antiquark pair, then the qq 
pair elastically scatters on the proton, and finally the qq 
pair recombines to form a virtual photon. This is shown 
schematically in Fig. The amplitude for the complete 
process is simply the product of the amplitudes of these 
three subprocess. We discuss each in turn. 

The amplitude for the incoming virtual photon (with 
polarization A = +,—,0) to fiuctuate into a quark- 
antiquark pair (of flavor / with helicities h,h = ±1/2) 
is given by the photon light cone wave- function {ijjyj)^'^ 
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which depends upon (Q^,z,r). Here z denotes the lon- 
gitudinal momentum fraction of the photon carried by 
the quark and r the transverse size of the qq pair. These 
wave functions are determined from light cone perturba- 
tion theory to leading order in the fcrmion electromag- 
netic charge |ir| . Similarly, the amplitude for the qq to 
recombine to form a virtual photon is (-0* ■ 

The elastic scattering of the qq pair with squared mo- 
mentum transfer — ~t is described by the elastic 
scattering amplitude, J^^(x, r, A). (The notation in this 
section follows S. Munier et al 5'].) Here Af^ is defined 
such that the elastic qq cross section is 



(1) 



It is convenient to work in coordinate space and define 
the S-matrix element at a particular impact parameter h 



5(6) = l-f ± / d2Ae^b-^A^f(a;,r,A) . (2) 



We denote amplitud.e for a vector meson to fluctuate into 
a qq pair by {^^vy^ ■ This vector-meson wave function 
will be discussed below. The amplitude for elastic diffrac- 
tive vector meson production with squared momentum 
transfer A^ = —t and transverse (T) or longitudinal po- 
larization (L) is 



Al,t(A) 



dz 
47r 



(9) 

Assuming that the S-matrix element is predominantly 

real we may substitute 2(1 — S{b)) with -grf^. Then, the 
elastic diffractive cross section is 



da 



L,T 



dt 



1 

16^ 



d^r 



dz 
47r 



d%{V\^),^^e 



-ib-A '^'^qq 
d^b 



(10) 

Note, it is the total dipole cross section at a particu- 
lar impact parameter which appears in this formula. S. 
Munier et al |1| extracted S{b) from the diffractive data 
using the formalism outlined above. 



This corresponds to ones intuitive notion of impact pa- 
rameter only when the dipole size is small compared to 
the size of the proton. The total cross section for the qq 
pair is found by taking the imaginary part of iA^i or in 
terms of the S-matrix element 



aqq{x,r) :^lmiAl^{x,r,0) ^ / 2(1 - Re S'(6)). (3) 



This motivates the definition of the qq differential cross 
section 



da, 



qq 



d^b 



2(1 -Rc S{b)). 



(4) 



The total 7*p cross section for transverse and longitu- 
dinal virtual photons is then given by the imaginary part 
of the forward scattering amplitude 

a0(x,g2) = jd?rj ^Y.^l\l)i^ <Jqq{x,r), (5) 



with 



aqq{x,r) ^ d b 



j2j, '^'^qq 



d% 



(6) 



III. MODEL DESCRIPTION 

In order to describe the data it is necessary to make 
specific assumptions about the dipole cross section and 
the photon and vector meson wave- functions. First, we 
give a simple model for the dipole cross section and con- 
trast this model with other models that have appeared 
previously. Then we define the gluon structure function 
xg{x,ix^) and the proton shape function T{b) which en- 
ters the dipole cross-section. Then we define and dis- 
cuss the properties of the vector meson wave-functions 
in some detail. Finally after the dipole cross-section and 
the wave- functions are specified, we determine the model 
parameters by fitting the inclusive total cross section and 
elastic diffractive vector meson data. 



A. The Dipole Cross Section 

The total cross section for a small qq dipole to pass 
through a dilute gluon cloud is proportional to the dipole 
area, the strong coupling constant, and the number of 
gluons in the cloud 



(7l7)t.T denotes the probability for a polarized virtual 
photon to fiuctuate into a qq pair with flavor / 

(7|7){ ^ Y.^r,j)f{^,j)f (7) 

hh 

(7|7)^^ ^ ^ E ir,j)f{r^,j)f- (8) 



hh 
A=±l 



^qq = -Y^'r'^ols{^J?)xg{x,^i'^) , 



(11) 



where xg{x,n^) is gluon density at some scale ji^ . Now 
imagine that the density of gluons in the target is not 
small. Divide the target into thin slices of thickness dz 
|46j |. The probability that a dipole at impact parameter 
b does not suffer an inelastic interaction passing through 
one slice of the proton is 



In the dipole picture elastic diffractive vector meson 
production appears in a similarly transparent way [ll| . 



P{b) = 1 - —r^as{n'')xg{x,^i'')p{b,z)dz . (12) 
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Here p{b, z) denotes density of the gluons within a proton 
and is normahzed to one 



(fibdz p{b, z) = 1 



(13) 



Exponentiating this result, the probabihty that the 
dipole suffers no inelastic interactions passing through 
the entire proton is |5'(6)P 



\S{b)f = exp 



N, 



■r'as{fi')xg{x,fi')Tib) .(14) 



Here T{b) is the thickness function, 

T{b) = / dzp{b,z) 



(15) 



To obtain the total cross section at a given impact pa- 
rameter we assume that the S-matrix element is predom- 
inantly real. Then the cross section at a given impact 
parameter b is given by 2(1 — ReS{b)) or 



dan 



d?b 



1 — exp 



2Nc 



r'^as{p^)xg{x,p^)T{b) 



(16) 

This expression is known as the Glauber-Mueller dipole 
cross section 'sj and can also be obtained within the 
McLerran-Venugopalan model JS] . It provides a sim- 
ple model for the interaction of a qq probe with a dense 
target. The function T{b) and xg{x^p?) are determined 
from the fits to the data. Figure |21 shows the dipole cross 
section as a function of the impact parameter b for var- 
ious dipole sizes at two distinct a;-values. The proper- 
ties of the dipole cross-section will be discussed below. 
Similar dipole cross sections were found in the model of 
Ref. 14]. The cross-section is evaluated for the largest 
dipole sizes only to illustrate saturation in the Glauber- 
Mueller approach. This dipole size is not contributing to 
the observed cross-section. 

The Glauber-Mueller dipole cross section is a simple 
model which respects unitarity. As discussed below, in 
the large region the model matches smoothly onto 
DGLAP since we will evolve xg{x,fj,^) with DGLAP. 
In the small x region, the model evolution can be im- 
proved. Indeed, evolution based upon the Balitsky- 
Kochegov equation [islflGjl has been extensively studied 
1(5, 7] and provides a reasonable description of the to- 
tal and total diffractive cross sections. Nevertheless, the 
Glauber-Mueller dipole cross section provides a simple 
model which exhibits saturation as a function of impact 
parameter in a reasonably generic way and is therefore a 
useful tool in exposing saturation effects in the data. 



B. The Proton Shape T{b) 

The impact parameter dependence of the dipole cross- 
section determines the properties of the distributions 
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FIG. 2: The dipole cross-section as a function of the impact 
parameter b for various dipole sizes evaluated at two x-values. 



in elastic diffractive processes. Here the variable t de- 
notes the squared four momentum transfer between the 
incoming and outgoing proton in the diffractive process, 
7*p Vp. The distribution is directly related to the 
Fourier transform of T{b) via Eas.1161 and 1101 Presently, 
the best data on diffractive <-distributions are available 
for vector-meso n p roduction. HERA experiments re- 
ported fvA IT^ ITU that all diffractive vector-meson t- 
distributions can be described by a single exponent B 
for \t\ < 1 GeV2 



dt 



cx exp(-B|i|) 



(17) 



Since the J/ip wave-function is better known than the 
p or (f) wave- functions, we will concentrate on J/^p pro- 
duction. The value of the coefficient B for the J/ip t- 
distribution is around 4 GeV~^ and is only weakly de- 
pendent and . The exponential form in Eq. El is 
the Fourier transform of 



Tcib) 



1 



2ttBg 



eM-bV^Ba) 



(18) 



We therefore tried this form for T{b). To account for the 
expected exponential behavior at large impact parame- 
ters we also considered the form 

TcYib)^ J dH'eM-i'^-^')V^BGY)Koib'/WGY) ■ 

(19) 
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FIG. 3: The differential cross-section for exclusive diffractive 
J/tl> production as a function of t for representative bins in 
W [13 • The solid (dashed) lines show the results of the IP 
saturation model assuming Tgy {Tg) for the proton shape. 



Due to the asymptotic behavior of Kq{x) ^ 
^/^^/2x exp{—x), TQY(h) falls off exponentially at large 
impact parameters. 

The parameters Bq^ Bqy and Wgy are only weakly 
dependent on the other parameters in the model and 
were determined iteratively. From the fit to the limited 
t range, \t\ < 0.7 GeV^, the following values were found: 
Bg = 4.25 GeV-2 for Tdb) and Bgy = 3.25 GeV-^ 
and Wgy = 0.66 GeV"! for TGvib)- Figure shows a 
representative subset of data [13 together with the re- 
sults of a fit to the data using both forms of T{b). The 
slope of the data fixes the parameters Bg and Bgy and 
Wgy ■ The normalization of the curves in Figure |31 is a 
determined by the underlying J/iJj light cone wave func- 
tion and is therefore model dependent. We will use three 
model wavefunctions which are discussed more fully in 
Sect, nil Fl and in the Appendix. The figure is for the 
CORNELL wavefunction. Unless explicitly stated, the 
profile function adopted throughout the rest of this work 
is TgyQ}) which for practical purposes is equivalent to 
Tcib). 



C. Gluon Structure Function and Dipole Evolution 

The dipole cross-section (Eq. [TB|) requires the gluon 
density xg{x, fj^ ) for all scales p?' . As in previous work , 
the scale fj? is related to the dipole size by 



C 



(20) 



The gluon density is evolved to the scale fj,^ with leading 
order DGLAP without quarks for the gluon density. The 
evolution of the gluon density is performed numerically 
for every dipole size r during the integration of Eq. 
The initial gluon distribution is taken at the scale /ig in 
the form 



xg{x,^il) 



A„ 



(1 



^5.6 



(21) 



The parameters of the model C, ^q, Ag, and Xg are de- 
termined from a fit to the DIS data. For light quarks the 
gluon distribution is evaluated at x — xbj — q2^^/2 ■ 
For the charm quark the gluon structure function is eval- 
uated at a: = (4m?J- Q'^)/ (Q'^ + W"^), where rric is the 
charm quark mass . This implies that the dipole cross 
section is flavor dependent when performing the sum over 
flavors in Eq. 13 

The present model has, in some sense, one less param- 
eter than the model of Ref . ^] . The parameter ctq which 
previously determined the asymptotic size of the dipole 
cross-section is now replaced by the parameter Bg (or 
Bgy and Wgy ) which is (are) fixed by the shape of the 
diffractive i-distribution. 

The gluon structure function obtained from the fit is 
shown in Fig. ^ for various dipole sizes. The dipole size 
determines the evolution scale /z^ through Eq. 1201 QCD 
evolution changes the rate of rise of the gluon structure 
function xg{x,fi^). The rate of rise is usually quantified 
by the exponent Ae//, 



A, 



eff 



d\ogixg{x,^j,'^)) 
d\og{l/x) 



(22) 



The exponent Ae/ / is a measure of the strength of gluon 
emission process, see e.g. |23 |. Figure shows a rapid 
increase of the exponent Ag/ / with decreasing x for small 
values of the dipole radius r. 



D. Comparison with the GBW model 

In the original GBW model the dipole cross section 
aqq{x,r) was independent of impact parameter and gave 
a good description of the inclusive total 7*p cross section. 
It also described some diffractive processes 'l',!2','2^ with- 
out providing any information about the i-dependence. 
In the GBW model the dipole cross-section aqq is given 
by 



GBW 

Oqq 



1-CXp ( 



(23) 
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FIG. 4: The gluon structure function for various dipole sizes. 
The dipole size determines the evolution scale /i^ . The dashed 
lines shows the GRV results 1211. 
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FIG. 5: The logarithmic rate of rise of the gluon structure 
function as a function x, X^ff — '"^^/p^^x/'^) ■ -^e// is a 
measure of the strength of the gluon emission process. 
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FIG. 6: The integrated dipole cross section as a function of 
the dipole size. The solid and dashed curves trace the dipole 
cross-section in the IP Saturation model and the saturation 
model of Ref. Q respectively. 



where Rq denotes the saturation radius and ctq is a con- 
stant. The saturation radius is x dependent 



1 



GeV2 \xq 



(24) 



The parameters ctq — 23 mb, Xgbw — 0.29 and xq — 
3 • 10~^ were determined from a fit to the data. The 
saturation radius Rg is analogous to the gluon distribu- 
tion and determines the growth of the total and diffrac- 
tive cross sections with decreasing x. For large size 
dipoles (jfg^^ saturates by approaching a constant value 
(To which is independent of the exponent Xqbw- 

Once the impact parameter is included saturation be- 
comes b dependent. For small values of 6, the cross- 
section grows rapidly with r until it reaches the satura- 
tion plateau seen in Fig. 13 d(jqq/<Ph = 2. The extension 
of the plateau bs can be estimated from the condition 
■^j^r'^as{iJ,'^)xg{x, fi'^)T{bs) = 1. Using the Gaussian 
form of T(6) and approximating xg{x,iJ,^) ^ {xo/x)^^ff 
we obtain 



b% 



2Bg log 



4NcBg 



(25) 



increases 



Thus, the integrated dipole cross-section (iqq 
logarithmically for large r and small x due to the growth 
of the plateau region. The logarithmic growth of the 
cross-section at small x is determined by Ag// and Bq- 
This behavior should be compared to the GBW cross 
section which becomes a constant uo that is independent 
of Xgbw- The differences are illustrated in Fig. El 
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E. Photon wave-function 



The longitudinal and transverse photon wave-function 
in the conventions of Ref. are given by 



Kojer) 
27r 



Koier) 



Koier) 



27r 



r 



(26) 



where 



Z(1-0)Q2 + ^ 



is the 



electromagnetic charge of quark flavor /, and drK^ier) = 
-eKi{er). 

For small dipoles (er <^ 1) the photon wave-function 
scales as i^i(er) ~ 1/er. For large dipoles (er 1), 
the wave-function falls off exponentially with the quark 
mass, Ki(er) ^ ■y/7r/2erexp(— er). 



F. Vector Meson Wave-Function 

To describe vector meson production, we developed 
several models of the J/ip light cone wave-function, 
{'4'v)'x^- These phenomenological wave- functions should 
satisfy several model independent constraints [2^1 l2^ . 
First, the wave-function must satisfy the normalization 
condition 



1 = 



(fkdz 
167r3 



hh 



ii'v)f\k,z) 



(27) 



Second, the decay width of the vector meson further con- 
strains the wave- function. Precise definitions and explicit 
formulas for the decay width and other constraints on the 
longitudinal and transverse wave functions are given in 
the Appendix. 

As an example, we model the J/ip wave-function as a 
Gaussian in momentum space in the meson rest frame 



ipNRip) = A^L,Texp(- 



-) 



(28) 



were p denotes the charmed quark 3-momentum. In DIS 
reactions the meson wave-functions are defined on the 
light-cone and not in the J/tp rest frame. Following 
Frankfurt and Strikman we boost the J/'f/' wave- 
function from the rest frame to the light-cone. The con- 
stants Nl^t and ujl,t were adjusted to satisfy the nor- 
malization and decay constraints. This wave-function is 
denoted Gaus-RF (Rest Frame). 

As noted in Ref. ,25] this procedure does not lead to the 
proper asymptotic form of the light cone wave function 



at small distances, (j) ~ -^(1 ~ z)- Therefore, we also 
considered a Gaussian wave-function defined directly on 
the light cone 



fc2 

4>L,T{k, z) = Nl z{l ~ z) exp(- — ) . 



(29) 



Once again, the norm and decay width constrain 
and lol- This wave- function is denoted Gaus-LC (Light 
Cone). 

J /ip wave-functions were also studied with potential 
models [23, |2^ and successfully describe the energy 
levels of cc bound states. Potential models were less 
successful at reproducing the meson decay width which 
is sensitive to the short distance behavior of the wave- 
function. This shortcoming may be less important since 
the short distance behavior of the wave-function is sup- 
pressed by the vanishing dipole cross section at small r. 
We determined the non-relativistic wave-function by inte- 
grating the radial Schrodinger equation with the Cornell 
potential [S^- As before we naively boosted this wave- 
function to the light-cone. Naturally, this wave-function 
is denoted CORNELL. 

In order to visualize the meson wave-function we define 
the meson photon overlap P7*y(?') as 



P-y-vir) = r / '^{i\V)t 



(30) 



which has the simple interpretation as the amplitude for 
diffractive vector meson production from the elastic scat- 
tering of a dipole of radius r. The overlap function is 
shown in Fig.dfor the transverse component of the three 
J /ip wave-functions considered in the paper. 



G. Dipole size 

To interpret the IP Saturation model physically, we de- 
termined the typical dipole size contributing to the total 



cross-section. The distribution of dipole sizes 



dr 
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is the integrand of Eq. ]E\ after the z integration is per- 
formed. Figure |S1 shows this distribution at representa- 
tive values of x and Q^. For each illustrated x and 
the arrow indicates the median radius of a dipole which 
scatters off the proton. Note that at smallest and 
X = 10^^ the majority of dipoles are smaller than r < 2 
GeV^^ and thus half of the dipoles are perturbative. 



IV. RESULTS FOR 7> PROCESSES 

A. Inclusive total cross-section 

The fit to the total 7*p cross-section data was per- 
formed in the range x < 0.01 and > 0.25 GeV^. The 
small X cut leads to an upper limit on Q^, < 200 
GeV^. Figure El shows the 7*p cross section as a function 
of at the various measured by the HERA exper- 
iments [s^ . Following the rise of the total cross section 
with W'^ at fixed Q^, we observe a striking change: for 
small (the top of Fig. ^ the cross section rises slowly 
with the center of mass energy W^, while for large 
(the bottom of Fig. ^ the cross section rises strongly. If 
we parameterize the cross section by a power law 

'^ir - (iy^)^-«2^) ^ (i/x)^-(Q^) , (31) 

this change in the rate of rise from low to high trans- 
lates into a dependence of the exponent Xtot ■ Mot (Q^) 
is illustrated in Fig. 1101 and summarizes one of the prin- 
ciple observations at HERA. At small Q^, cr^.p increases 
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FIG. 7: A comparison of the meson photon overlap function 
(Eq. I3UII for the three different J/ip wave-functions adopted 
in this work. 



FIG. 8: The distribution of the dipole sizes contributing to 
the inclusive total cross-section at various virtualities and 
some representative x values. The vertical arrows indicate 
the median radius of a dipole which scatters off the proton. 
The horizontal arrows indicate saturation region r > rs(b,x) 
at b = 0. See Table [TTI 

with as in all hadronic reactions, Xtot ~ 0.08 [sj. At 
at large the exponent Xtot increases substantially. 

In the IP Saturation model the increase of the ex- 
ponent Xtot with is a consequence of the DGLAP 
evolution. When the dipole size becomes smaller than 
~ ^ DGLAP evolution becomes increasingly effective 
and the growth of the gluon distribution at small x be- 
comes increasingly rapid. Mathematically speaking, Ae// 
increases for small size dipoles as seen in Fig. O Aver- 
aging over dipole sizes with the photon wave function 
translates Ae//(r) into Xtot{Q^)- The large increase in 
Xtot with indicates a transition from a predominance 
of large to a predominance of small dipoles when evalu- 
ating the cross-section. The contribution of large dipoles 
at low is determined by the light quark mass, mu,d,s- 

The parameters Ag, Xg and /ig were found by fitting 
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FIG. 9: The 7*p cross section as a function of at various 
. The values of are shown on the left side of the plot 
along with the scale factor applied to the data for a better 
visibility. The full (dashed) line shows a fit by the IP Satura- 
tion Model with mu,d,s = 50 MeV (m„,d,s = 300 MeV). The 
fits are performed for x < 0.01. 



mu,d,3 (MeV) 


A, 


A. 


2 

Mo 


X^/Ndf 


300 


1.67 


0.226 


6.85 


702/160 


150 


2.34 


0.043 


1.32 


199/160 


100 


2.85 


-0.04 


0.96 


153/160 


50 


3.47 


-0.12 


0.77 


137/160 



TABLE I: Table of parameters determined in the fit to the 
total 7*p cross-section for different assumptions on the light 
quark mass mu,d,s- 



the total 7*p cross-section for various light quark masses. 
The parameter C is strongly correlated to the parameters 
Ag and Ag and does not affect the fit quality. Therefore, 
we set C = 4 by fiat. The mass of the charm quark was 
rUc = 1.3 GeV as required by elastic diffractive J/ip pro- 
duction discussed below. The results of the fits to the 
ZEUS data are collected in Tabled Because of the good 
agreement between the HI and ZEUS data sets the com- 
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FIG. 10: The exponent \tot as a function of . 
parametrizes the total cross section at fixed , crj 
(l/s:)^*°*. 



bined Hl-ZEUS fit returns nearly identical parameters. 
The fits show a clear preference for the light quark mass 
of mu^d^s = 50 MeV, which is taken as a reference value. 
Note, the quality of the fit reported here is even better 
than the fits of Ref. 0. 

As seen in Table U when light quark mass is increased, 
the fit subsequently increases Ag. DGLAP evolution can 
only increase the observed rate of rise of the cross section 
with W'^, Xtot > Ag. Therefore, as seen in Fig. O for 
larger masses of the light quarks the cross section rises 
too quickly with at low Q^. It is certainly known 
that the strange quark mass is heavier than the u and d 
quarks. However, the influence of the strange quark on 
the overall fit is relatively small since its contribution to 
F2 is suppressed relative to the u, d contribution by the 
ratio of charges squared, es/(e„ -I- e^) = 1/5. Therefore, 
we have neglected the difference between the s and u, d 
quark masses. Including a heavier strange quark would 
make the fit only marginally worse. However, the clear 
dependence of the fit quality on the mass of the light 
quark is a property of the model. It is due to the interplay 
of DGLAP evolution and the properties of the photon 
wave function. It may indicate the shortcomings of the 
wave function at large distances. 

In contrast to the GBW model ^ and the model of 
ReL 2], saturation is not responsible for the change in 
Xtot with once the impact parameter is included. 
This is because the total dipole cross section continues to 
grow within the saturation region although the differen- 
tial dipole cross-section has attained its asymptotic value 
^ = 2. Therefore, for the total cross-section saturation 



9 



effects are smaller than in the original GBW model. How- 
ever, for diffractive ^-distributions saturation will appear 
more clearly. 



B. Charm Production 
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FIG. 12; The ratio Fl /-F2 as a function of x at various photon 
virtuahties. 
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FIG. 11: A comparison of the measured F2 |32| to the results 
of the model. 

The dipole model makes a direct prediction for the in- 
clusive charm contribution to F2 . In the dipole approach 
the charm quark distribution is calculated from the gluon 
distribution. Figure ITTI shows a comparison between the 
measured and predicted values of the charm structure 
function F2. The results depend weakly on the charm 
mass. The good agreement with data for both ZEUS 
and HI experiments 'sS, 's^l confirms the consistency of 
the model and supports the dipole picture. 



C. Determination of Fl 

As a byproduct of our investigation we determined the 
longitudinal structure function F^. A measurement of 
Fl would be a critical test of our understanding of the 
gluon distribution. F^ is zero in the naive parton model 
and is non-zero only after gluon radiation is included. In 
dipole models at small x, both Fl and F2 are governed 
by the gluon distribution and therefore they behave simi- 



"3d 



10 



IP Sat 

MRST 


r=0.5 GeV"' 

r=l GeV"' 

r=2 GeV"' 

1-3 GeV"' 

r=S GeV"' 











10 



10 



10 



10 
X 



FIG. 13: A comparison of the gluon structure function of the 
IP Saturation model with the results of the MRST fit ^ for 
various dipole sizes. The dipole size determines the evolution 
scale fj,^ . 
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larly. Figure El shows the ratio of to F2 as a function 
of X for various photon virtuaHties. The ratio is around 
18 % and gently varies with x. In the standard QCD 
fits Fr / Fi may show a stronger variation with x at low 
[33 ■ This is because the shape of the gluon struc- 
ture function seen in Fig. 1131 is very different from the 
dipole model presented here. The difference in the gluon 
structure function between the MRST fit [s^l and the 
dipole model fit is due to the sea quark contribution. In 
the MRST fit the initial sea quark contribution is inde- 
pendent of the gluon distribution and subsequently influ- 
ences the gluon evolution. In the dipole fit the sea quark 
distribution is always derived from the gluon distribution. 
It is important to measure Fl in order to determine the 
gluon density precisely and to clarify the relative roles of 
gluons and sea quarks. 

D. J/ijj production 
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FIG. 14: A comparison of the measured total elastic diffrac- 
tive J/ip cross-sections [l^ with the results of the model. 

The cross sections for elastic diffractive J/ip produc- 
tion were calculated using Eq. EH with the dipole cross 
section determined by fitting aj*p. Figure FPil compares 
the computed total J/ip cross-section with HI measure- 
ments [13 for various photon virtualities at = 90 GeV. 
The absolute magnitude of the cross-sections is very sen- 
sitive to the charm quark mass. The charm quark mass 
was adjusted to nic = 1.3 GeV in order to match the 
experimental cross section at = 0. A small shift of 
50 MeV in the charm mass changes the cross-section by 
around 15%. 

The predicted total J/tp cross-sections for the Gaus- 



FIG. 15: A comparison of the measured total elastic diffrac- 
tive J/ip cross-sections |0 with the results of the model. 



RF and Gaus-LC wave-functions are very similar in the 
range considered. Therefore, we only show the predic- 
tions of the Gaus-RF wave- function in what follows. The 
CORNELL wave-function, which is more concentrated at 
the origin leads to a prediction which is in significantly 
better agreement with data at large Q^, see Fig. E] In 
contrast to Ref. , the scale fi^ = C/r'^ + is fixed 
with C = 4. Within the range C = 1 ~ 10 the results 
were almost independent of the parameter C after refit- 
ting total 7*p cross-section. 

Figure 1151 compares model predictions with the mea- 
sured J/ip cross-section as a function of energy W at 
= 0. In this case, both J/tp wave- functions are in 
good agreement with the data. 

Note, there is a ~ 15% discrepancy in the normaliza- 
tion of the J/tp cross-sections between ZEUS and HI. 
This discrepancy was absorbed in the charm quark mass; 
we take rric = 1.25 GeV when comparing with ZEUS data 
and rric — 1.3 GeV when comparing with HI. 



E. t distributions for J/tp production 

Figure 1161 compares the J/tp diffractive cross-section 
as a function of t for various to model predictions. 
In the measured region, both J/tp wave- functions lead 
to very similar results although at large t they start to 
differ significantly. In the large t region the exponential 
behavior of the differential cross section is considerably 
distorted and the IP Saturation model predicts diffractive 
dips. Figure [T7I illustrates the emergence of the dip by 
displaying the Fourier transform of the differential dipole 
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FIG. 16: A comparison of the measured difFractive J/xp cross- 
sections fl9| with the results of the IP Saturation model for 
two different J/-)/) wave-functions. Also shown are the results 
of the non-saturated model evaluated with CORNELL wave- 
function. 



cross section for various dipole sizes 



d<7„ 



dt 



(32) 



For small r, the h dependence of the dipole cross section 
is essentially Gaussian and its Fourier transform is a pure 
exponential in t. For large r, due to saturation, the shape 
of the dipole cross section becomes more box-like, as seen 
in Fig. 121 The two dimensional Fourier transform of a box 
is proportional to the Bessel function Ji which leads to 
dips in the t distribution. Although the distortions to 
the proton shape seen in Fig. |2]are small when displayed 
in coordinate space, they appear vividly after taking the 
Fourier transform and examining the t-distribution. For 
dipoles with r > 3 GeV~-^ the i-distribution is substan- 
tially changed. 



F. Saturation versus Non-Saturation 

In order to isolate the effects of saturation we compared 
the results of the IP Saturation model with a model which 
is not unitarized. This model, denoted IP Non-Sat, re- 



FIG. 17: The Fourier transform of the differential dipole 
cross section for various dipole sizes. 



places the dipole cross-section of Eq. ^2 with the first 
term after expanding of the exponential 



don 



d% Nc 



(33) 



The proton shape T{h) was assumed to be the same as 
in the saturated case. The fit to the total cross sec- 
tion was repeated and a non-saturated gluon structure 
function was determined. The optimal fit was obtained 
for a higher quark mass than in the saturated case, 
iTLu^d^s = 150 MeV. This is expected since the dipole 
cross-section aqq grows faster with r than in the saturated 
case. The fit yields Ag = 3.5, \g = -0.13, = 0.80 
and /Ndf = 0.8. Naturally, the saturated and non- 
saturated dipole cross sections are very similar since they 
both fit the inclusive 7*p cross-sections very well. Conse- 
quently, all observable cross-section which only are sen- 
sitive to the dipole cross-section integrated over impact 
parameter are indistinguishable. This was verified explic- 
itly for the total diffractive J /ip production. In particu- 
lar, this agreement can be seen in Fig.^jfor the Xtot{Q^) 
distribution. However, the gluon distributions are not the 
same as illustrated in Fig. ^| which compares the satu- 
rated and non-saturated models. For small dipole sizes 
the distributions are similar but they differ significantly 
as r becomes large. 
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The t-distribution provides additional insight into the 
dynamics of saturation in DIS. Figure [TCI compares pre- 
dictions for difFractive J/V' production in the saturated 
and non-saturated modes . As discussed above saturation 
distorts the impact parameter distribution and leads un- 
avoidably to difFractive dips in the t-distribution. With- 
out saturation, as seen in Fig. 1161 the exponential de- 
crease continues to large t. Diffractive dips in the t dis- 
tribution were clearly observed in elastic hadronic reac- 
tions [sgIS . In DIS processes at HERA the measured range 
of the ^-distribution is too small for a dip to appear. In 
addition the HERA experiments have always measured 
elastic diffractive events with some small admixture of 
inelastic diffractive events. Even if inelastic events are a 
small fraction of the total, these events will dominate the 
i-distribution at large t since the inelastic t-distribution 
falls slowly ^E^. 

In the small t range, saturation slightly increases the 
slope parameter B of diffractive vector meson production 
at large W. Figure El compares the slope parameter B 
determined by the ZEUS experiment with model pre- 
dictions. The observed increase of the slope B with W 
(called shrinkage) is slightly larger than the model predic- 
tion. This could be an indication that the shape T{b) is 
somewhat dependent on x due to QCD evolution. Since 
considerably more data is available than was used in the 
present analysis we plan further study of the observed 
shrinkage. 

Concluding we note that the i-distribution provides 
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FIG. 18: A comparison of the gluon structure function in 
the IP Saturation and IP Non-Saturation models at various 
dipole sizes. The dipole size determines the evolution scale 




FIG. 19: The diffractive slope B measured by the ZEUS 
experiment as a function of W at = 0. The solid line 
shows the prediction of the IP Saturation model. 

precise information about the interplay of saturation and 
evolution phenomena. Therefore, it is crucially impor- 
tant to measure elastic vector meson production with a 
small systematic uncertainty over a wide range in t. 

V. SATURATION SCALE g| 

In the limit of high gluon density called the Color Glass 
Condensate (CGC)_4J, the interaction probability in DIS 
becomes large. The probability of no interaction is mea- 
sured by the square of the S-matrix element, S^{b). Nor- 
mally, in high energy reactions the interaction probabil- 
ity is small, S"^ ~ 1. However, in the CGC limit this 
probability approaches unity, S'^ ^ 0. Figure [201 shows 
S'^{b = 0) as a function of x for various dipole sizes. 
For small size dipoles the gluon structure function grows 
strongly with decreasing x and the S-matrix diminishes 
quickly. For large size dipoles, the gluon density density 
does not grow and the S-matrix is small even at relatively 
large x. 

To delineate the saturation region in the {x^r'^) plane 
we first define a density profile D 

Dib, X, r') = ^as{fiHr^))xg{x, fi^{r^))T{b) . (34) 

Then, the saturation radius rs is defined as the dipole 
size for which proton consists of one interaction length 

5^=exp(-:^^^)=e-. (35) 
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FIG. 20: The square of the S-matrix as a function of x for 
various dipole sizes. 



FIG. 21: The saturation scale as a function of x. The dashed 
curve shows the critical line (S'^ — e^^) in the Non-Saturated 
model at 6 = 0. 



This defines an implicit equation for r| which can be 
solved numerically. Now the saturation scale (5| is de- 
fined as the density profile D evaluated at the saturation 
radius 



Qlib,x) = Dib,x,rl)^—. 



(36) 



In the GBW model S"^ = exp(-rV2i?2) and we find 
QI = 1/^0 accord with the original definition of the 
saturation scale. The saturation region is then defined 
by the condition that r > rs{b,x). 

Figure 12] shows the saturation scale as a function 
of l/x in the central impact parameter region. We ob- 
serve that the saturation scale is lower than the GBW 
value. This is primarily because the charm quark con- 
tribution to F2 was neglected when the saturation scale 
was evaluated. Including the charm quark reduces the 
gluon density and the saturation scale. 

We also calculated the S-Matrix element and deter- 
mined the corresponding effective saturation scale for 
the non-saturated model. As seen in Fig. inn the sat- 
uration scale is very similar for IP Saturated and IP 
Non-Saturated models. This is because the correspond- 
ing gluon structure functions are very similar for small 
size dipoles, as seen in Fig. ^| Therefore, our estimates 
of the saturation scale are only weakly model dependent. 

Table lU records the saturation radii and the percent- 
age of the total 7*p cross-section in the saturation region, 
r > rs- At low (Q^ = 0.4 GeV^) the cross-section is 
dominated by dipoles which are larger than the satura- 
tion radius. As seen in Fig. a fraction of these dipoles 




FIG. 22: The saturation exponent As as a function of x and 
impact parameter. The dashed curve shows the exponent 
along the critical line (5*^ — e~^) in the Non-Saturated model 
at b=0. 
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may be considered perturbative since they are smaller 
than r < 2GeV~"'^. In the large region saturation is 
of secondary importance. The region k, A GeV^ is 
particularly interesting since a significant fraction of the 
dipoles are perturbative and are approximately the size 
of the saturation radius. BFKL and saturation dynamics 
are expected to influence observable cross sections in this 
kinematic domain. 

It is also interesting to deflne the saturation exponent 
As 



d\og{Q%) 
dlog(l/x)- 



(37) 



Figure shows the saturation exponent As as a func- 
tion of 1/x for various impact parameters. The sat- 
uration exponent exhibits a sizable x dependence. It 
varies between As = 0.15 at cc = 10~^ and As = 0.35 
at a; = 10~^ which should be contrasted with the con- 
stant GBW value, Xgbw — 0-29. An analysis of BFKL 
evolution near the saturation boundary found that the 
exponent As is nearly constant as a function of log(l/a;) 
in the HERA kinematic window 37]. 
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FIG. 23: The transverse density ATA(h) for several light 
nuclei compared to the proton transverse profile, Tp{b). 



VI. NUCLEAR DEPENDENCE 

Now that impact parameter distribution of the proton 
is taken into account, the dipole picture provides a frame- 
work for modeling nuclear phenomena in DIS. The distri- 
bution of nucleons in the nucleus pa{^ is parametrized 
according to the Woods-Saxon distribution [38| 



PA{r) 



N 



exp 



(r-B.A) 



(38) 



where A denotes the number of nucleons in the nucleus, 
5 = 0.54 fm, Ra = {l.l2bn) A^/^ - (0.86 fm) ^-^/^ and 
N is adjusted to normalize the distribution to one 



(frpAir) = 1 . 



(39) 



b (GeV-1) 





1 


2 


(GeV^) X 


rs % 


rs % 


rs % 


0.4 10^^ 


1.1 72 


1.3 68 


2.0 52 


0.4 10"'' 


1.6 70 


1.9 65 


2.7 51 


4.0 10"'' 


1.6 33 


1.9 25 


2.7 17 


4.0 10"^ 


2.2 44 


2.4 28 


3.2 20 


40. 10"^ 


2.2 13 


2.4 12 


3.2 9 


40. 10'^ 


2.7 21 


2.9 18 


3.5 16 



TABLE II: Table of saturation radii and of the percentages of 
the inclusive total 7*p cross-section in the saturation region, 
r > rs- 



The transverse distribution is defined as 

f + oc 

TA{h) = 



dzpAW\P + z2 

-OO 



(40) 



Figure [53 compares the proton shape Tpib) with the 
transverse density TA{b) for several light nuclei. 

For a first estimate of nuclear effects we assume that 
the gluons are distributed uniformly throughout the nu- 
cleus. This assumption is realized by replacing the proton 
shape in Eq.^jby the integrated Woods-Saxon distribu- 
tion Ta(6) scaled by the number of nucleons, 



da„ 



<Pb 



1 — exp 



2Nr_ 



,{l?)xg{x,t?)ATA{b) 



(41) 

This naive assumption, called the smooth nucleus as- 
sumption below, cannot be true since Fig. 1231 shows that 
the center of the proton is on average more dense than 
the carbon nucleus. Nevertheless it is instructive to eval- 
uate /AF^ for this case and to compare it with NMC 
measurements of shadowing at small-x [s^, see Fig. 1241 
The figure shows that the smooth nucleus assumption un- 
derestimates nuclear shadowing considerably as can be 
anticipated by comparing the densities in Fig. 1231 

The proton profile Tp{h) indicates that the gluons are 
distributed within individual nucleons with a radius of 
w 0.5 fm. The nucleons are then distributed within the 
nucleus according to a Woods-Saxon distribution. There- 
fore the gluons form a lumpy distribution within the nu- 
cleus. To account for this correlation among the gluons 
we generalize the model in the following way. The prob- 
ability of a given configuration of nucleons within the 
nucleus {(bi, zi), (b2, ^2), • • • , (bA, za)} is 

ntip(bi,z,)- (42) 
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FIG. 24: Shadowing in DIS: a comparison of model predictions (open points) to data (full points) with the smooth and 
lumpy nucleus assumptions. The NMC measurements |39| were performed at = 0.74, 1.20, 1.70 GeV^ for Carbon and at 
= 0.60, 0.94, 1.40 GeV^ for Calcium (in order of increasing x values). 

For a given configuration of nucleons the probability that the dipole scatters at impact parameter b is 
I 



daqq 



= 2 



{(bi,2:i),(b2,2:2),--- ,(bA,ZA)} 



1 - exp [ -^r'^as{n^)xgix, y?) ^ Tp{h - b;) 



i=l 



Novif we may average over all possible configurations of nucleons to determine the cross section 



(Pb 



(fbi dzi - ■ ■ d^bA dzA -P{(bi,zi),(b2,22),- - ,(bA,ZA)} 



daqq 
IPb 



{(bl,2i),(b2,Z2),---(bA,2^)} 



Noticing that the integrals factorize we have 

2[1 - (1 - I)^] 



d% 



I = / dH'TA{h + h) 



1 — exp — 



2Nc 



r 



(43) 

(44) 

(45) 
(46) 



To simplify this further, we notice that the integral 
over b is different from zero only over the size of the 
proton which is small compared to the nucleus even for 
the deuteron. Thus TA(b + b') « TA(b). Using the 
definition of the total cross section of the qq pair on the 
proton aqq (Eq. Il)land[]i)| we have 



daqq_ ^ 
d^b 



1 - (1 - 



(47) 



Only if A is large and r is small do we recover the smooth 
nucleus formula of Eq. 1411 Similar considerations were 
discussed in the context of setting up initial conditions 
for classical QCD simulations of nucleus-nucleus colli- 
sions Figure 121 shows that the lumpy nucleus model 
slightly overestimates the observed shadowing. However, 
good overall agreement with data clearly indicates that 
the IP Saturation model can be used to study of nuclear 
effects. 
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To estimate the saturation scale for heavy nuclei we 
follow the same procedure as for proton. We set the 
lumpy nucleus S-matrix element squared, jS'p — (1 — 
TA(b-o) ^p_'^2A Qq^oi This is an implicit equation 

for the saturation radius, rs- Once rs is determined 
numerically the saturation scale is determined from the 
relation, Ql — 2/r|. 

Figure 1251 shows the saturation scale for various nu- 
clei calculated in this fashion. It is interesting to ob- 
serve that the saturation scale is not increasing like A^^^. 
Rather, for Au the nuclear enhancement factor is only 
a factor of w 2 at x = 10~^ and « 3 at a; = 10"''. 
When comparing the saturation scale to other estimates 
it is important to remember that this is the satura- 
tion scale for quarks. The saturation scale for gluons 
is {Ql)g = ^iQl)g « 2.25 (Q2)^ [111. Thus for RHIC 
conditions with x « 10~^, {Ql)g is « 1.3 GeV^ which is 
slightly smaller than was estimated previously by analyz- 
ing RHIC multiplicity data IH lil IH Q . 
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FIG. 25: The quark saturation scale {Ql)q at 6 = as a 
function of x for various nuclei . The gluon saturation scale 
is (Q')9 = ^(Q?), ^2.25(Q2),. 



VII. SUMMARY 

We developed an Impact Parameter dipole saturation 
model which describes HERA the total 7*p cross sec- 
tion and differential diffractive J/ip with high precision. 
The new feature of the model is an impact parameter 
dependent dipole cross-section which respects unitarity 
constraints. The impact parameter distribution is deter- 



mined by measured J/tp ^-distributions to a large extent. 
For the functional form of the dipole cross-section we 
adopt the Glauber-Mueller formula which can also be 
obtained within the McLerran-Venugopalan model ^3 . 

The parameters of the dipole cross-section are deter- 
mined from a fit to HERA data. The fit depends on the 
specific form of photon and the J/ip wave function. At 
short distances, the photon wave-function is well known 
from QED. For the J/ip wave- function we proposed sev- 
eral ansatzes and verified that the dipole cross-section 
was almost independent of these alternatives. 

QCD evolution effects are taken into account by evolv- 
ing the gluon structure function with the lowest order 
DGLAP equations. DGLAP evolution is essential in con- 
trast to the GBW model. 

In the dipole approach only the gluon structure func- 
tion is modeled. All sea-quark contributions to F2 are 
obtained from the photon wave-function and the subse- 
quent dipole-gluon scattering. The validity of this ap- 
proximation is corroborated by the excellent agreement 
between the predicted and observed open charm contri- 
bution to F2. 

The impact parameter dependence of the dipole cross 
section determines the properties of the J/tp total and 
diffractive cross-sections. The saturation ansatz for the 
dipole cross section predicts diffractive dips in the t distri- 
bution which are not necessarily predicted by an ansatz 
without saturation. For J/tp production the diffractive 
dips appear in a region of t which is not covered by data 
but which is accessible to future measurements. In ad- 
dition, the saturation ansatz predicts a slight increase of 
the diffractive slope B, which reflects a small increase in 
the proton size at small x. The change of the slope B with 
energy, called shrinkage, provides information about the 
interplay of saturation and evolution phenomena. Since 
shrinkage was observed for several vector-meson species 
and since considerably more data exists than was used in 
the present paper, further study of diffractive vector me- 
son production and their t-distributions is warranted. We 
also point out that additional measurements of vector- 
meson production in a extended i-range would provide 
crucial information about the gluon density. 

With an accurate determination of the S-matrix ele- 
ment we evaluated the saturation scale - the scale 
at which multiple collisions becomes important. The Ql 
determined in the model is approximately half of the 
GBW value and is only weakly dependent on the spe- 
cific form of the dipole cross-section. 

In the low region up to ~ 4 GeV^, saturation dy- 
namics contributes significantly to the observed cross- 
sections. The saturated state is at least partially per- 
turbative since it is frequently initiated by small dipoles. 
Above ~ 10 GeV^ saturation plays a secondary role. 
Dipoles are very small and although they are strong gluon 
emitters the gluon phase-space density remains small. 

In order to isolate the effects of saturation, we also 
studied a dipole model which is not unitarized . This IP 
Non-Saturated model also obtained a good description of 
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the data. However, even with the non-saturated model 
the fit returned large gluon densities and unitarity cor- 
rections became important in the same region of x and 
. Thus, the saturation ansatz is the more consistent 
approach. 

The extension of the model to nuclei is fairly straight- 
forward. We evaluated the nuclear cross-section with two 
different ansatzes for the distribution of gluons within 
the nucleus and compared the results to data on nuclear 
shadowing at low x. We found that the lumpy nucleus 
ansatz, in which gluons are distributed around the nucle- 
ons in a similar manner as around a free proton, leads to 
a good description of the data. We then used this ansatz 
to evaluate the saturation scale in heavy nuclei and found 
that Ql increases more slowly then A^/'^. Nevertheless, 
for Gold and Uranium is enhanced by factor of '-^ 2 
when compared to the proton. For RHIC conditions the 
gluon saturation scale was found to be « 1.3 GeV^, 
which is somewhat smaller than was estimated previously 
from RHIC multiplicity data jU] . 

The success of the model suggests that the saturated 
state known as the Color Glass Condensate plays an im- 
portant role at HERA. Saturation is particularly impor- 
tant in the small x and low region where the 7*p total 
inclusive cross-section exhibits the same universal rate of 
growth as all known hadronic cross-sections. Revealing 
the Color Glass Condensate experimentally would repre- 
sent a fundamental advance towards an understanding of 
the high energy limit of QCD. 
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APPENDIX A: THE J/t/; WAVE-FUNCTION 

The discussion here closely follows 26] and uses some 
of the concepts of . The light cone wave- function for 
J can be written as an expansion in the Fock space, 

\V) = ^Pcc |cc) + i^qqg \qqg) + .... (Al) 

Below we neglect higher Fock components in the wave- 
function and treat the meson as a predominantly a cc 
pair. Restoring indices, the J/tp state is given by 



(Pkdz , ~ 



Sab 



167r3 



(V^y)f (k, Z)^ \c'^{k, Z, h) c'i-k, 1 - Z, h)) 



(A2) 



Here A = 0,±1 denotes the polarization of the vector 
meson, h,h = zt denote the helicities of the c-quark and 
anti-c-quark respectively, and a, b are the color indices of 
the c-quark and anti-c-quark. Here and below we strictly 
adhere to the light-cone conventions of Ref . . 

The J/ijj wave-function is modeled after the photon 
wave-function. The longitudinally polarized vector me- 



son (A — 0) is 



^h+h- 



h+ 



V2 



4>L{k,z), (A3) 



while the two transverse polarizations (A — ±1) are given 
by 



J 



{^v)X\{k,z) 
{^vt\{k,z) 



+ — e^"^'{z5,,+j,_ - (1 - 2:)5^_,- + ) + 5,,+,- + 



e-^^-iil - z)S,,^j,_ - zS,,_j,^) + 6^,_j,_ 

TTlr 



<i>T{k,z) 
4>T{k,z), 



(A4) 



where fce"^'" — kx + iky. Once the wave-function in mo- the wave- function in coordinate space 
mentum space is specified then Fourier transforms return 

(^y)f (r,z) = J ^e*-(^v)f (fc,z) . (A5) 
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Note under Fourier transform ±ke^^'^'' (f){k,z) becomes function is then 
=pi e^"'''' z). The overlap with the photon wave- 



(7|f)l^E(^;)o'(V'v)5 

hh 



2z(l-z)Q^<p, 



hh 
A 



hh 
A=±l 



2tt 



r 



1 



Two inputs are used to constrain the wave-function 
(j)L,T- The first input is the normahzation of the wave- 
function, (v^' \V^) = {2Tif2p+5{p'+ - p+)S^{v' - 
\i)5^ ^ ■ This leads to the requirement that 



1 = 



cPkdz 
167r3 



hh 



cPkdz 
167r3 



(A7) 



(AS) 



for 



and 



dPkdz 
167r3 



(z2 + (l-z)2) +1 



(A9) 



rric dr 



•T ^ eKi (er) 



2n 



Pt rUc 



Kojer) 
2tt 



(A6) 



for (f)T ■ Note that the formula (B19b) in [23 is incorrect 
by a factor of two for the transverse wave- function. The 
second input is the leptonic decay width of the vector 
meson. The decay width of the meson is given by 



where fv is given by 



My 



(AlO) 



(All) 



In this formula = (e^,e ,ex,£y) is the polarization 
vector of the meson and e is the charge of the positron. 



and 



1 / \/2(0, 0, 1, ±i) for transversely polarized mesons 
e(; = 0,0) for the longitudinally polar- 

ized mesons. Evaluating this matrix element in light cone 
perturbation with the conventions of I3| gives 



fv 
fv 




4'T{k,z) 
z{l-z) 



(A12) 



where Cc — +| (in contrast to Cc — \/4Tra^ec) is the 
charge in units of the positron charge. 

To go further, we must specify the form of the wave- 
function. We have studied three model wave- functions, 
GAUS-RF, CORNELL, and GAUS-LC. Two of the mod- 
els (GAUS-RF and CORNELL) are motivated by the ex- 
pected long distance behavior of the wave-function. One 
of models (GAUS-LC) is motivated by the expected short 
distance behavior of the wave- function. 

For large distances, and for truly heavy vector mesons, 
the non-relativistic wave-function should give a good 
description of the light cone wave- function. To make 
the correspondence between the light cone and non- 



relativistic wave-functions functions we repeat the dis- 
cussion in Ref. 25]. In the non relativistic limit, terms 

proportional to (^:^^ can be dropped and the the nor- 
malization condition Eg. lATl is identical for the longitudi- 
nal and transverse wave- functions. Simple kinematic re- 
lations give the invariant mass of the the quark antiquark 

2 I .2 

pair My = (pi + P2)^ = T(i-z) • This relation assumes 
that the quark and antiquark are on shell, which is true 

for non-relativistic quarks to order - In the rest 

frame of the meson (where the quarks have momentum p 
and — p) the invariant mass is {pi -f ^2)^ = 4(p^ + m^). 
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Comparing these terms gives the identification 
fc2 + (2z - \fml 



V 



4z(l - z) 



d kfdz . 



(A13) 
(A14) 



Requiring probabihty conservation per unit phase space 



(27r)3 



-ipNRip) 



2 d^kfdz 



167r3 



(t)L.T{kt, z) 



(A15) 



leads to the identification 



4>L,T{kt 



L,T 



\Az^{l-zf^ 



fpNR P 



Ikl + (2z- l)2m2 
4z(l - z) 



(A16) 



Here ipNR.ip) is the non-relativistic wave-function in mo- 
mentum space 

i^NR{x) = j ^e^P-^^NRip) . (A17) 

For the two long distance models (GAUS-RF and 
CORNELL) 4'nr{p) is specified. Then with the form 
of 4'L,T specified, the non-relativistic motivation for the 
light cone wave- function is abandoned and Nl,t are con- 
stants adjusted to satisfy the normalization and decay 
conditions (Eqs. IA7I and |A12|I including the terms pro- 
portional to {kt/rric) ■ Terms proportional to (kt/rric) 
mock up relativistic effects and constitute « 25% of 
the total elastic diffractive cross section. For the short 
distance model (GAUS-LC) 4'l,t is specified directly. 
For the numerical results quoted below we take My — 
3.096 GeV and = 1.3 GeV. 
GAUS-RF . In this model we take 

^NR.{p) = A^L.rexp J— ^ ' (^18) 

where ujl.t and A^l,t are adjusted to satisfy Eq. IA12I 
for the decay width and Eq. IA7I for the norm. These 
parameters are Nl = 1.74 GeV"^/^, ujl = 0.66 GeV and 
Nt = 1.86 GeV-3/2 and ujt = 0.61 GeV. 



CORNELL . In this model we take the form of tpNii{p) 
from the wave-function of a non-relativistic potential 
model of J/ip. Using a shooting method solve 
the radial Schrodinger equation for the Cornell poten- 
tial to find iJjNRjx). The parameters of the potential are 
taken from Ref. [23. Note, the constituent quark mass 
of the potential model {nic = 1.8 GeV) differs from the 
current quark mass {rric — 1.3 GeV) in Eqs. IA3IIA4IIA16I 
for the the light cone wave-function. In this model the 
decay condition Eq. IA12l is not satisfied. 

GAUS-LC . We now turn to the short distance model 
GAUS-LC. At short distances the light cone vector meson 
wave-function is proportional to cx z(l — z) This 
motivates the following ansatz for the wave-function 

(t)L.T{r,z) ^ Nl,tz{1- z) eM —) (A19) 

Nl^t and ujl.t are adjusted to reproduce the normaliza- 
tion condition and the decay width. For J/ip we find: 
Nl = 1.21 GeV-\ lul = 0.80 GeV and Nt = 1.63 
GeV-i and ujt = 0.58 GeV. 

With the wave-function (pL^xik, z) specified the wave- 
function (j)L,T(j', z) and its derivative dr(j){r, z). are found 
taking Fourier transforms. The final overlap between the 
photon and vector meson wave-function is given by Eq. 
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